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A method of synthesizing a general nth order phase-locked loop is pre- 
sented. In contrast to conventional phase-locked loops, the circuitry is 
digital rather than analog. The general circuit consists of an assembly of 
logic blocks (gates and storage elements) which, when driven by external 
clock signals, exhibits phase-locked loop properties. These properties, along 
with high stability and the absence of adjustments, make the digital phase- 
locked loop ideally suited for use in large systems which use monolithic 
integrated circuits for microminiaturization. Analysis and synthesis tech- 
niques make use of Z-transform methods in achieving the desired frequency 
response as the realization of an nth order difference equation. A general 
technique is developed and two specific cases, n = 1 and n = 2, are con- 
sidered in detail. Analytic results relating to the phase-locked loop's static 
and dynamic performance are derived and found to correlate well ivith 
laboratory results for actual circuits. 

I. INTRODUCTION 

A new phase-locked loop (PLL) with interesting properties has 
been developed for potential application in large multiple data set 
installations which provide low speed serial data communications 
for time-shared computers. An objective for such data set arrange- 
ments is to minimize cost per channel by putting the major part of 
the required circuitry into a common section where it may be shared 
by all channels. This objective is achieved by using a digital PLL 
as an FM demodulator with low cost logic circuits located in the 
channel units and clocks with their associated driving amplifiers 
located in the common circuits. PLL's which use analog circuits have 
received considerable attention and analysis and synthesis methods 
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are available. 1 - 2 However, the circuits covered here are digital, and 
the approach is similar to that of digital (or sampled-data) filters. 8 - 4 
By using a digital PLL, no low-pass filter or voltage-controlled oscil- 
lator, generally associated with the feedback loop of conventional 
PLL's is required.* This property, along with high stability and the 
absence of adjustments makes the digital PLL ideal for microminia- 
turization using monolithic integrated circuits. 

This paper presents snythesis procedures for an nth order digital 
PLL. The PLL realized by such a procedure possesses a response 
which obeys a linear nth order difference equation. Analysis is per- 
formed using Z transform methods commonly encountered in sampled- 
data control systems. 5 The technique is that of establishing a mapping 
between the s plane and the z plane so that a correspondence between 
the coefficients of the controlling nth order difference equation and 
the desired s plane poles may be established. An iterative circuit 
is presented so that once the coefficients are determined, the nth 
order loop may be realized. 

The remainder of the paper is devoted to the synthesis, realization, 
and analysis of two specific examples, n = 1 and n = 2. Both systems 
are analyzed to determine static and dynamic performance, and the 
results of data transmission tests are given. The out-of-lock behavior 
as well as internal noise resulting from jitter is characterized for the 
first order PLL. The second order PLL is representative of higher 
order systems and the analysis is easily extended. The "capture 
phenomenon" associated with underdamped systems is encountered. 
In both of the examples considered, experimental results are found 
to correlate well with theory. 

n. THE DIGITAL PHASE-LOCKED LOOP 

As a prelude to the synthesis procedure we show that the loop 
response of the PLL can be expressed as an nth order difference 
equation. Figure 1 is a block diagram of the loop. Among its basic 
components are an exclusive or comparator which develops an output 
gating function dependent upon the phase relation of its inputs, and 
a transmission gate acting on several clock signals /«, gi to provide 
inputs to register circuitry. When the loop is locked, a shift circuit 
periodically transfers the contents of the (i — l)th register to the 
ith register; the period is one half that of the input signal. The shift 

♦Although no internal loop filter is needed, a low-pass filter is required to 
recover the demodulated baseband signal. 
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Table I — Symbol 

a, Real part of z-plane pole, z,-. 

0i Imaginary part of z-plane pole, z,-. 

7, 5 Parameters pertaining to stability analysis of second order PLL. 

A Normalized total input frequency deviation. 

A/ Total input frequency deviation. 

A//, Lock range of PLL. 

Av Voltage step. 

Aw.v Peak-to-peak noise voltage. 

A^sig Peak-to-peak signal voltage. 

*u, eo', «i Time errors resulting from phase discontinuity when switching 

between clocks. 

m Real and imaginary parts of Butterworth characteristic with carrier 

frequency input. 

£(*) Duration of kth positive (logical "1") level of feedback (flip-flop) 

signal. 

p(k) Interval between A;th and (A; + l)th zero crossing of input signal. 

94 Real part of s-plane pole, s,-. 

r{k) Duration of kta positive (logical "1") level at exclusive- or output. 

u ( Radial cutoff frequency of low-pass filter, 

w, Imaginary part of s-plane pole, s,-. 

C Counting capacity of JV-stage counter. 

E Input voltage. 

Ef Feedback voltage. 

E PLL output voltage. 

/ Frequency of input signal. 

fa, ft Discrete input frequency. 

fc Cutoff frequency of low-pass filter. 

fm Average input (carrier) frequency. 

fr Rest frequency of PLL. 

/«, ft End point of frequency lock range. 

/., 0. Clock frequencies (» = 1, 2, • • - , n); also denotes input to registers. 

Fi, Gi Normalized clock frequencies (i = 1, 2, • • • , n). 

H(s) Transfer function (Butterworth). 

i . <-«*: 

hji k, I Counting indices. 

M Threshold of final register. 

n Order of system; number of registers. 

N Number of counter stages per register. 

N m Minimum N for stable operation. 

3j jth s-plane pole. 

/ Time. 

Td Period of baseband data signal. 

T\, T«, T 3 , T t Signal durations for stability analysis. 

v(k) Average output voltage during fcth interval 

v(l) Continuous time function. 

v 2 Average output voltage for parasitic mode. 

z Argument of z-transform. 

Z[v(k)], V(z) Z-transform of v(k). 

(and reset of the first register) is controlled by the nth register which 
provides an output pulse after the Mth clock pulse is counted. A flip- 
flop converts the output pulse train to a square wave and provides 
an input to the comparator. 

The difference equation describing operation of this circuit is de- 
veloped with the aid of waveforms shown in Fig. 2, which gives the 
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Fig. 1 — Block diagram of general digital PLL. 



steady-state terminal signals for the comparator. The output signal 
is the gating function. Assume that clocks g lt g 2 , 0s, • • • > 9n are 
enabled during the "0" level, and all other clocks j x , f 2 , . . . , U are en- 
abled during the "1" level. With the first register initially cleared, its 
count at the conclusion of the (fc + l)th period of gating is g x [p{k) - 
T (fc)] + /n-(fc + 1). During each successive period, this count is 
shifted into the ith register and augmented by a count of gi[p(k + 
i _ i) - T (fc + % - 1)] + f iT {k + i) for i = 2, 3, . . . , n. The count 
propagates through the n registers where, upon reaching the number 
M at the nth register, the count is reinitiated. Although n periods are 
required for a complete count cycle, the process may be thought of as 
the interleaving of cycles initiated p seconds apart. 
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III. GENERAL FORM OF THE DIFFERENCE EQUATION 

It can be seen from Fig. 2 that the count of the nth register is the 
sum of pulses counted during the intervals p{k) — r(k), P {k + 1) — 
r(k + 1), . . . , P (k + n - 1) - r{k + 11 - I), t(/c + 1), T (fc + 2), 
. . . , r(k + n). Summing the counts for the sources fi, fa, . . . , f n , 
0i j 92 , • • • , Q n find equating the sum to M gives 

OMQ ~ r(k)] + g 2 [p(k + 1) ~ r(k + 1)] + • • • 

+ g n [p(k + n - 1) - r(k + n - 1)] 
+ UrQz + 1) + • • • + f n r(k + n) = iV, (la) 
or rewriting, 

/„T(fc + n) + (/„_, - 0„)r(A; + n - 1) + • • • - ^r(fc) 

= M - [g nP (k + n - 1) + <7„-,p(A- + n - 2) + ■ ■ • 
+ g 2P (k + 1) + y lP (ft)]. (lb) 
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Fig. 2 — Input, feedback, and output waveforms for PLL which is assumed to 
be locked with constant input frequency. Gated clock signals are shown 
symbolically. 
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This is the system difference equation relating the response r{k + i) 
to an excitation P {k + i). In its general form this equation, when 
Z-transformed, is analogous to a Laplace transformed system equation 
in which polynomials in s multiplying the response and excitation 
functions result in poles and zeros, respectively. Since there is a 
one-to-one correspondence between the s and z planes with regard to 
poles and zeros, it is expected that by properly choosing coefficients, 
equation lb may be synthesized to provide a desired frequency re- 
sponse (high-pass, low-pass, bandpass, and the like) possessing spe- 
cified critical frequencies. As we show, it is the objective of this paper 
to exploit the low-pass properties of equation lb. To achieve maximum 
high frequency attenuation, the coefficients g ( ,i = 2,3, . . . ,n will be 
set equal to zero thereby locating all z-plane zeros at infinity. 

/„r(fc + n) + f n . lT (k + n - 1) + • • • - gAk) = M - g lP (k). (lc) 

Equation lc is normalized by letting Fj = fj/f„, j = 1, . . . , n — 1 
and Gx = fifi//„- Also, notice that the cycle-by-cycle average voltage 
of the r(k + j) interval expressed as a fraction of the maximum 
possible voltage is given by 

* + *-$$* < 2) 

Incorporating these substitutions gives 

v(k + n) + F n - X v{k + n - 1) + • • • + F x v{k + 1) - G,»(fc) 

= 2MF(k) - Q u (3) 

where 

F(k) = 



2f nP (k) 

If it is assumed that P (k) changes very little with k and that p(k) 
is small with respect to the system response time (which it is) , the 
v(k) can be represented as a sampled continuous function of time, 
v(t), letting t = k P . That is, v(t) is a function whose value at the 
fcth zero crossing of the input signal is v (fc) . 

This type of equation is best solved using z transform methods 
assuming that the input frequency is constant. In particular, the 
following transform pairs are noted : 

Z[v(t)) = Z[v(k)] m V(z) = fXjpV. ( 4a ) 
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Z[AexpW] = _^_. (4b) 

Z[v(k + j)] = 2 '[7( 2 )] - £ sV^ 1 - 0- (4c) 

i=i 

Accordingly, equation 3 is transformed as follows: 

+ F n _ l2 n -'[F( 2 )] - **„_, £ sy« - 1 - /) + • • • 

1-1 

+ F l2 [7(2) - »(0)] - G,7« = (2il/F - Cwfrfrj • (5) 
Combining terms, 

= (2W - Gof-^-r) + Z E ^'Ki - /). (6) 

where 

F n = 1. 

IV. GENERAL SYNTHESIS PROCEDURE 

Before developing a synthesis technique for the digital PLL, it is of 
interest to review a specific application, that of FM demodulation. A 
"lock range" may be defined for the PLL whereby steady state input 
signals having constant frequencies lying within this range will cause a 
steady-state output, v(k + j) = v(k + ; + 1), all ;, such that this output 
is linearly related to the input frequency / = l/2p. This is the relation 
required for demodulation. 

For dynamic behavior one may consider a binary baseband signal in 
which each of the two states is assigned a discrete frequency /„, j b within 
the lock range. If the baseband signal switches randomly between states 
with a maximum rate \/T d a new phenomenon is introduced. In this 
context the PLL may be regarded as a low-pass filter which should 
possess a bandwidth, u e , equal to or greater than ir/T d . By proper choice 
of the coefficients in equation 6 the desired filter shaping may be a- 
chieved. It is expected that the low pass filter characteristics will be a 
function of the input frequencies /„, j b . The resulting complication is 
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conveniently eliminated by making use of the narrowband approxima- 
tion /„, f b Psi i„ where / m = $(/«. + U)- This is not unrealistic, because the 
PLL was developed for just such a narrowband system.* With this in 
mind, we now give a synthesis procedure. In using Z transform techni- 
ques it is assumed that the input frequency is approximately constant 
so that samples are taken at equal time intervals. 

The coefficient of V(z) in equation 6 is the "characteristic poly- 
nomial" of the system and using it, the desired low pass filtering 
properties of the loop can be synthesized. For example, assume that 
a transfer function with poles in the s plane at Si , s 2 , . . . , s„ is to be 
synthesized. Its characteristic polynominal is given by 

n <•-«,). en 

A conformal mapping between the s plane and z plane is given by 
the transformation 

z = exp(ps), (8) 

by which the pole s,- = cr, + «o,- is mapped into 

z t = exp [GO (<r, + ia { )] 

= [exp (po-,)] [cos (co,p) + i sin («,-p)]. (9) 

Since the complex s plane poles occur in conjugate pairs, and sin (w,p) is 
an odd function, the corresponding z plane poles also occur in conjugate 
pairs, and the desired characteristic equation is transformed into 

jl (z- z,) = z\+ A^z"- 1 + • • • + A v z + A, , (10) 

i-i 

where each of the coefficients are real. Equating coefficients in equa- 
tion 6 to those in equation 10 gives the required clock frequencies 
0i, fi, ...,/«• A negative value implies an associated register which 
counts "down" while a positive value suggests a register which 
counts "up." 

V. FHIST ORDER DIGITAL PLL 

5.1 Static and Dynamic Behavior 

The difference equation for an n = 1 PLL is easily written from 
equation 3: 

* Full duplex transmission with } a = 1070 Hz, /> = 1270 Hz in one band, and 
} a = 2025 Hz, U — 2225 Hz in the other band. 
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v(k + 1) - Giv(k) = 2MF - G„ (11) 

in which M = 2 N ~ X and N is the number of counter stages in the 
register. For a steady state condition to exist, t>(fc + 1) = v{k) so that 



v(k) = 



2MF - G t 2Mf - g x 



(12) 



1 - &\ /. - 0, 

Since ^ v(k) ^ 1 the end points of the lock range are given by 

/ I ,(»-o - U = gi/2M; / | .(„., - /„ - U/2M. (13) 
And so the lock range is 

1 



Af L = \iu - ii\ = ojj? I & ~~ h I- 



(14) 



The static response is diagrammed in Fig. 3. 

The dynamic response to a step change in frequency is given by 
equation 6 for n = 1. 



[2 - G l ]V(z) - [2MF - ©J ^jj + z»(0). 
A partial fraction expansion yields 



(15) 
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Fig. 3 — Average output voltage versus input frequency for first order PLL. 
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Assume that for time t < the normalized input frequency is F a . 
At t = the input steps to F b . Then from equation 12: 

m - W±^- (17) 

Substitution into equation 16 gives the z transform of the response 
for t ^ 

The inverse transform is easily found with the aid of equation 4a. 

... 2MF b - G x 2M(F h - F.) {f1 , k 
v( ® = 1-G, ~ 1-G l {Gl) 

m 2Mf b -g x _ 2M(j b - f a ) faY for k ^ m 

/i - Q\ 7i - Qi W 

Assuming that v(k) is a continuous function of time, we let v(k) = 
v(t) and t = k P = k/2f h . This gives 

„ (/) = m* - * _ 2J/(/> - /■) exp [_ 2 / fc /][ln C/i/^i)]. (20) 

The resulting time constant is 

T - (21) 

1 2/ 6 In (h/ 9l ) ' K } 

and thus the half-power bandwidth is 

/„ = -/6ln(/ 1 /^ 1 ). (22) 

7T 

Notice the dependence of filter shaping upon input f b . It follows that 
there is a somewhat different time constant for input frequency 
changes from j b to f a . 

A simplified first order digital PLL has been built as shown in Fig. 
4. It is a special case of the general system of Fig. 1. The FM input 
was composed of the discrete frequencies j a = 1070 Hz and f b = 1270 
Hz. Clock frequencies g x and ji were chosen as (2M) (970) Hz and 
(2M) (1370) Hz, respectively, to provide a lock range of twice the 
total input frequency deviation. M was chosen to be 128 (see Section 
5.3). The above values gave a time constant, T, of 1.13 ms. Inter- 
changing f a and f b resulted in an increased time constant of 1.34 ms. 
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Fig. 4 — Block diagram of simplified first order PLL FM discriminator. 

Oscilloscope photographs showing the two step responses of the actual 
loop are shown in Fig. 5. 

5.2 Out-of-Lock Oscillation 

If there is no input to the PLL, it will run at a rest frequency of 



7r \/, "*" gj " M gi +U 



(23) 



Notice, however, that one half cycle will be at Q\/2M and the next 
at/i/2M. 

If an input to the PLL is present, but its frequency lies outside 
the lock range, that is, if / < /j./ 2M or / > g 1 /2M } the output voltage 
v{t) will oscillate between and 1. Figure 6 shows waveforms for a 
loop which is out of lock. For time to the left of the dotted line, the 
loop is attempting to lock and the frequency of the feedback signal, 
E f , is approaching that of the out-of-lock input E. Within this region, 
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Fig 5 — Output voltage response of fust order PLL to input step frequency 
of (a) 1070 to 1270 Hz and (b) 1270 to 1070 Hz. Horizontal scale: 0.5 ms per cm. 

equation 16 applies and is rewritten in the time domain as 

For time to the right of the dotted line in Fig. 6, a new difference 
equation applies. It is written as 

gi[p(j) - r(j + 1)] + fiTtf) = M (25) 

so that a derivation similar to that previously used results in 

Equations 24 and 26 are represented graphically in Fig. 7. M , gi , and 
/i take on the values of the previous example and an out-of-lock 
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Fig. 6 — Input, feedback, and output waveforms of nonlocked first order PLL. 

input frequency / = 1400 Hz is chosen as the input. The positive 
slope segment to the left corresponds to equation 24 with an initial 
voltage v(0) assumed to be zero. When v(k) reaches unity, the next 
segment is governed by v(j) with an initial condition v'(0) = 1. 
When v{j) reaches zero, the response is again given by v(k) and with 




TIME, t 



Fig. 7 — Single cycle of output average voltage oscillation for nonlocked first 
order PLL. 
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Fig. 8 — Actual output voltage of experimental first order PLL showing out- 
of-lock oscillation. Horizontal scale: 1 ms per cm. 

initial condition adjusted accordingly. The actual waveform for the 
experimental PLL using this example is shown in Fig. 8. 

5.3 Internal Noise From Phase Discontinuity 

The clock signals are unsynchronized, and thus a random phase 
discontinuity results at the time of gating. This results in an internal 
noise voltage which exhibits itself as a fluctuation in the output voltage. 
This effect can be thought of as a quantization noise of the phase. A 
portion of Fig. 2 is redrawn in Fig. 9 to show the time errors c , t' Q , and e, 
at the instant of gating caused by differences in clock phases. It is 
assumed that the frequency of the ripple is low enough so as not to be 
attenuated by the PLL filtering property. With this in mind, difference 
equation 11 may be amended to account for the phase error. 

[r(* + 1) - ei(k + 1)]A 

+ [p(k) - r(k) - e (k) - e' (k)) gi + 2 = M. (27) 

For a constant input frequency and very slowly varying r(k), r(k + 1) = 
r(k) so that 

r(k) [/,-<?,] = M - g lP + faiQc + 1) 

+ gxMk) + «*(*>)] - 2, (28) 



or 



[»(*)] l/i ~ 9i\ - 2M/ - * + 2// l€l (/t + 1) 

+ 2fo,[€o(fe) + <(k)) - 4/. (29) 
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&v N is defined as the difference in voltage for maximum and minimum 
(zero) phase errors. Thus 



Since 
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e,(fc + 1) | M + g*Mk) | m „ + e' (k) | max )]. (30) 



ei(fc + 1) In... = T and e (fc) | mBS = e' (k) | m „ = — , 
7i <7i 



then 
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(31) 
where / is approximated by the carrier /„, for the narrowband case. 
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Fig. 9 — Waveforms for locked first order PLL showing how time jitter is 
caused by phase differences in clock signals. 
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It may be seen that the peak-to-peak output for an incoming FM 
wave with phase continuity in the clock signals is given by equation 
12 as 



Ay„ iB = 



™1« ~ fc _ 2M 1>- ~ g' I = 2M A/ , (32) 

/i — 0i U — d I /i — 01 ' 

where A/ is the total frequency deviation \f a — /• |. Thus the mini- 
mum peak-to-peak voltage signal-to-noise ratio within the loop is 

^ = 2MM = 2 „£. {33 ) 

' Avat 6/ m 6/ m 

Continuing the previous example with A/ — 200Hz and / m = 1170 Hz, 
it is desired that the internal noise be 20 dB below the signal so that 
2 N ^ 361. This suggests a nine stage counter (2 9 = 512); however, it 
should be pointed out that a worst case of phase jitter has been assumed 
and the average phase jitter is less than this. Laboratory experiments 
have shown satisfactory results with an eight stage counter thus giving a 
minimum signal-to-noise ratio of 17.2 dB. 

It is apparent from equation 33 that the internal noise may be made 
as low as desired by choosing a sufficiently large value of N. Thus a 
direct relationship exists between equipment cost (number of counters) 
and performance (jitter distortion). 

VI. SECOND ORDER DIGITAL PLL 

The location of the desired s-plane poles must be specified and mapped 
into the 2-plane. The following example is concerned with synthesizing 
the familiar Butterworth response, but the procedure is certainly 
applicable to other filter classes. 

A second order PLL is to be synthesized so that its response is that of 
a Butterworth low-pass filter with cutoff at oi e . The 2n = 4 poles of 
H(s)H(—s) lie on a circle of radius u e and subtend equal arcs such that 
the filter's characteristic equation is 

s 2 + (2)*(M + co 2 = 0, (34) 

with poles at 



«,.,-fijH±0. (35) 



The corresponding poles in the z plane are 

z li2 = exp (s li2 p) = [exp (-a)] [cos (a) ± i sin (a)] 
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where 

Digressing for a moment, it is of academic interest to map the entire 
Butterworth circle into the z plane. This is easily done by letting s be a 
circle of radius w c , that is, s = a, exp (?'<£), ^ <£ ^ 2ir. Substituting 
2 = exp (sp) and approximating / by j m gives 

z = jexp ||k- cos (</>) Jjj cos (2j- sin (0)J + i sin y~f- sin (0)^ 

(37) 

The resulting cardioid-like shape is shown in Fig. 10 with /„//« as a 
parameter. The angle </> is also shown so that portions of the Butter- 
worth circle may be conveniently translated into the z plane. 

Returning to the specific example, the characteristic polynominal is 
written using equation 36 as 

(z - z,)(z ~ *•) = 2 2 - [2 exp (-a)] [cos (a)]z + exp (-2a). (38) 

The resulting coefficients are equated to the respective coefficients in 
equation 6, namely z 2 + FiZ—Gi so that: 

F, = /,//, = -2[exp(-a)][cos«] (39) 

and 

Gi = gi/U = -exp (-2a). 

A third equation is arrived at by fixing the steady-state voltage v(k) 
for a specific input frequency, /. The choice is arbitrary, and since the 
system's input spectrum is symmetrical about the carrier frequency, 
f m , it is reasonable to fix the corresponding output voltage at 0.5. The 
equation is obtained from equation 6 by noting that steady-state im- 
plies periodicity so that v(k) = v{k+j). Thus: 

v(k)[l + F x - G 1 ] = 2M *= - (?, . (40) 

Solving equations 39 and 40 gives 

_ 2il//„ exp (a) 
n cos (a) - sinh (a) ' K } 

_ 4il// m COS (a) ( . 

71 cos (a) - sinh (a) ' K } 
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Fig. 10 — Mapping of the s-plane Butterworth circle into the 3-plane. 



2Mj m exp (-«) 
cos (a) — sinh (a) 



(43) 



As a practical example, consider the demodulation of a narrowband 
FSK wave whose spectrum is centered about a carrier of 2125 Hz. The 
signaling rate is limited to ^ 300 baud so that a cutoff frequency of 250 
Hz should prove adequate. As previously noted, the choice of M is de- 
pendent upon the maximum baseband jitter (quantizing noise) which 
can be tolerated. This must be weighed against the added circuitry and 
higher clock frequencies imposed by large values of M. A convenient 
choice is M = 128. Thus 

a = Wu = 0362 
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so that 

g, = 598 kHz, 

/, = 1504 kHz, 
and 

U = -1012 kHz. 

As previously mentioned, the negative value for / 2 requires a "down 
counter" for the second register. This example is continued in sub- 
sequent sections. 

6.1 Step Response of the System 

Once the coefficients g x , /i , . . . , /„ are determined, the characteristic 
equation is uniquely specified and the system's time response to a 
step in input frequency is easily found. The step response of a second 
order system is now derived in the interest of providing further in- 
sight into the characteristics of the loop. 

Equation 6 is written for n = 2 as : 
[z 2 + F lZ - O l ]V(z) 

= (2MF - G,)(^-y) + v(l)z + v(0)z 2 + zF lV (0). (44) 

Assume that the input frequency steps from F a to F at t = (as before, 
capitalization denotes normalization so that F, = /,// 2 , G t = f/,// 2 ). 
Thus for t > we have 



m - wf - q ,)[ + n _ Gl)(z n y] 



+ v(0) 



y + f x z - 



j + [M0) + KD][ 22 + ^_J - (45) 



The initial conditions v(l) and v(0) are determined from the differ- 
ence equation 

v(k + 2) + F x v(k + 1) - GMk) = IMF - G x . (46) 

Assume that the system is in steady state for t ^ so that v(k) = 
v(k - i), k ^ 0. For k = -2 

w(0)[l + F, - fly = 2MF a - G u (47) 
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so that 

Ml 2MFa ~ gl 

«(o) - 1 + Fi _ 0i - 

Now let fc = — 1 : 

«(1) = 2i¥F - G, - [F, - fly »(0), 
which may be written as 

y (l) = „(o) + 2MA, 

where 

A = [F - F a ). 
Substituting equations 50 and 51 into 45 and simplifying leads to 

2MA 



(48) 

(49) 
(50) 
(51) 



T 2J1/A 11" 

" Li + f, - gJL- 2 



+ F t - GUI/ + F lZ - Gi 

r gMAg, i r z 1 

" li + f, - gay? + F iZ - gA 



(52) 



The following substitutions are made so that V(z) may be easily 
transformed. Let 

Fi = -2[cos/3] exp (-a); Gi = -exp (-2a). (53) 

Substituting and rearranging gives 

2MA 



2MA 1 [" z 2 - z(cos |8) exp (-a) 1 

1 + F, - G: JU 2 - 2z(cos j8) exp (-a) + exp (-2a) J 



4MAG, 



- F, (-4G, - F,)' 



_\1 + F, - G, 

z(sin |8) exp (—a) 



* _z 2 - 2z(cos 0) exp (-a) + exp (-2a) _ 
This is transformed into 

v(k) = [v(0) + x + ™ A _ Gl _ u W 



(54) 
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— exp (— ak) 



2MA 



_\l + F, - G 



cos (0k) 



+ ( l + M Jfl Gi ~ ^i)(-4Gi - Fd~*Bm 08fc)j , 



(55) 



where u(k) is the unit step. The time response is found by again assum- 
ing that v(k) is a continuous function of time v(t), k « 2jt. Substitution 
gives the familiar response of an underdamped second order system: 

2MA 



v(t) = \v 



<0) + 



1 + F l - G l _ 



u{t) 



- exp (-2a//) ( - + £*_ G ) cos (2/3//) 



(56) 

The experimentally determined step response of the actual n = 2 
PLL closely approximates the Butterworth response and is shown in 
Fig. 11. 

6.2 Lock Range for the Second Order System 

As was the case for the n = 1 PLL, the loop is said to be in "lock" 
if, for a steady input frequency, the output v (k) is constant. Although 
there are many frequency ranges for which the loop exhibits this 
property, there is a "fundamental" lock range over which v(k) varies 
linearly with frequency between its defined limits, and the output 
frequency equals the input frequency. In steady state, v(k+2) = 
v(k+l) = v(k) so that from equation 40: 

2MF - G, _- 2Mf - g x 
f. + fi - 0. 



v{k) = 



l + F i ^G 1 

The end points of the lock range are 



if g v(k) ^ 1. (57) 



J |»(fc)-o — Ju — 



and 



The lock range is given by 



2M 

/. + u 

2M 



(58) 



(59) 



AU= \1«-h\ = 



2M 



0i - /i - U 



(60) 
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f i =2225 HZ — 



fj =2025 HZ 




— f; =2025 HZ 



YoutW 



•fi =2225HZ 



Fig. 11 — Output voltage response of second order PLL to input step fre- 
quency of (a) 2225 to 2025 Hz and (b) 2025 to 2225 Hz. Horizontal scale: 1 ms 
per cm. 

Substituting the values calculated previously gives a lock range of 
414 Hz which is approximately twice the bandwidth required. The 
useful bandwidth is reduced still further as will be shown. 

6.3 Bound on Cutoff Frequency 

For a step in frequency /<,-/» we have a corresponding step in volt- 
age, 



Ay = 



_ 2M/ a - gl _ 2Mf b - gl = 2M(ja ~ U) 



(61) 



U + fi - 0i U + 1\ - 9i U + 1i - yi 
Av is defined to have a maximum value of unity so 

/. + /i - ft fc 2M(/ a - U) (62) 
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and 



1 + F t - G x £ ^ (/ a - U) for U < 0. 

)2 



(63) 



If f a and /„ are symmetrically distributed about the carrier frequency, 
f c , we can let v = % at /i = /„, so that from equation 40 



U = 



4M1. 



(64) 



(65) 



1 + F, + <7, 

Substituting equation 64 into 63 gives 

1 + F l - ft £ (1 + F t + G r o(^p i )- 
Substituting equation 53 for F a and Gi gives 
1 - 2[exp (-a)] cos /3 + exp (—2a) 

- (S^) [1 " 2(C ° S ^ GXP ( ~ a) ~ GXP ( ~ 2a)] - (66) 

The constants « and /? are related to the filter shaping desired and are 
a function of the input frequency for t > 0. For the case of a Butter- 
worth response (n = 2) 



a = = 



(2)* \f 



» 



For a narrowband system the input frequency may be approximated 
by the carrier f m . Thus 



1 — 2(cos n) exp (— n) -f- exp (— 2/i) 

£ [1 - 2(cos M ) exp (- M ) - exp (-2 M )][|£J , 



(67) 



where 



This may be written as 






exp (n) - 



cos n ^ 



JbT 2 



7m g 



exp (-/*) 



(68) 
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Thus for a given f a , f b and f m , the cutoff frequency f is bounded. For 
example let f a and f b be 2025 Hz and 2225 Hz, respectively, and be 
symmetrically distributed about the carrier /„, . Then 

cos ^ exp(M) + (0.91)exp(- M ). (flQ) 

Solving this transcendental equation gives 

U ^ 0.19 (70) 

and 

j e ^ 182 Hz. 

An upper bound on f c is found by requiring that the end points of 
the lock range be positive frequencies. The lower edge of the lock 
range is required to be greater than zero so that from equations 59, 
41 and 42 

^ t/ 8 > or 2 cos M - exp (m) ^ 0. (71) 

This is satisfied if /x ^ 0.54 which gives positive value for equation 71 
as well as for the upper edge of the lock range g t /2M. Thus 0.19 ^ n ^ 
0.54 so that for f m = 2125 Hz 

182 ^ f, ^ 542 Hz. (72) 

The bounds on f /f m are loose in the sense that it is assumed that the 
full lock range is available to input frequencies. This is not the case 
for underdamped systems, where a hysteresis effect known as "cap- 
ture" reduces the effective lock range available, and so the bound 
might be tightened accordingly. 

6.4 The Capture Phenomonen • • 

When the PLL is "out of lock" its output voltage oscillates be- 
tween 1 and 0. The "pull-in" frequencies are those frequencies furthest 
removed from the carrier for which the PLL will ultimately lock. In 
general, the PLL will exhibit a hysteresis so that the pull-in range 
will be smaller than the lock range; also, the upper and lower pull-in 
frequencies will generally not be symmetrically distributed about the 
carrier /,„. In general, solution of the capture range (for conventional 
phase-locked-loops) results in a nonlinear integrodifferential equa- 
tion. Solution for conventional PLL's requires phase plane techniques 
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and is documented in Ref. 2. An analysis to determine capture range 
for the digital PLL has not been performed. 

6.5 Stability of the Second Order System Within the Lock Range 

Thus far there has been no restriction placed on N, the number of 
counters required per shift register. For N less than some critical 
value N m it is possible for stable modes of operation, different from that 
of Fig. 2, to exist. These modes do not exploit the PLL to its fullest 
advantage and thus in previous derivations it has been assumed that 
N ^ N m . 

A possible parasitic mode is shown in Fig. 12. The input frequency 
is constant, as is the frequency of the flip-flop signal, E f , but the 
waveform of E f is no longer square. This operation is caused by the 
limited length of the registers. Assume that the lengths (number of 
counter stages) for the n = 2 system are equal, and of value N. Each 
time the feedback flip-flop (see Fig. 1) is triggered, register 1 is reset 
to zero. During the period Tx+T 2 (see Fig. 12) this register is ad- 
vanced to a count of [T 1 g 1 +T 2 fi] modulo C where C = 2 N . This num- 
ber is shifted into the second register where one of two conditions may 
exist: 

T l0l + TV, ^ M, 



-p— 



Ef 



T, *-T,^ <- 



3 I 
Fig. 12 — Waveforms for parasitic locked mode of second order PLL 
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or 

The first equation results in normal operation, analyzed in previous 
sections. During T 4 the count is reduced by a negative / 2 so that 

[TxQt + r a /i + T<f a ] mod C = M. 
The second equation requires that the second register count down to 
zero (which is congruent to C) and then further reduce the count by 
C -M. Thus 

[T l g l + T 3 U + T4 % \ mod C - M - C. 
These equations may be combined by introducing the parameter 8 = 
— 1, depending upon the mode of operation: 

T l9i + IT,/, + ZV, - M - 8C. (73) 

Similarly, an equation similar in form to equation 73 covers the in- 
tervals T 2 , T 3 , and T 4 : 

T agi + T<U + ZVs = M -yC, (74) 

where y = 0, 1. 
The period, T 2 , may be explicitly solved for by noting that 

r, + r 4 - t* + t 3 = p. (75) 

Thus 



_ T, g, 

1,2 ~ P " /i + / 2 - * 



, (-8C + M)U - (~yC + M)(f 2 - gl ) . 



(76) 



Under normal conditions 8 = y = so that 



2Mj - gi 

V2 ~ U + U-g,' 



which is identical to equation 57. The parasitic mode of operation is 
given by8 = — 1, y = 1 so that equation 76 may be rewritten as 

U + U - gx /i + Tt - ffi 

Graphical examination (see Fig. 13) of equations 57 and 77 shows that 
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2M 
NPUT FREQUENCY, f 



Fig. 13 — Comparison of parasitic and normal mode voltage versus frequency 
characteristic for second order PLL. 

two modes of operation may exist within the lock range. To insure that 
this will not happen, C must be chosen so that the frequency range for 
which ^ v 2 ^ 1 lies outside the normal lock range. Thus at the upper 
lock range edge, / = g 1 /2M, it is required that v 2 > 1 so that from 
equation 77, 



2C 



Si. 

2M 



/. + gi - U 



> l. 



Hence 



C_^ /. + g. - U 



M 



> 



and thus the minimum number of counters required is 



N m = 1 + integer value 



(78) 



VII. PERFORMANCE 



In a laboratory test, an FM digital signal with mark and space fre- 
quencies of 1270 and 1070 Hz was fed into an n = 1 PLL which had 



* That is, the integer value of 23 is 2. 
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an 8-stage counter and clock frequencies for which g x = 2M(970)Hz 
and ji = 2M (1370)Hz. The output of the exclusive-OR was fed to an 
amplifier which clipped the signal to precise levels and provided a con- 
stant output impedance. The signal then went to a Butterworth low- 
pass filter with n = 3 and cutoff of 200 Hz. The filtered signal was sliced 
to give a digital output which could be compared with the original 
digital input. The FM modulator used in these tests was from a Bell 
System Data Set 103E1. The eye pattern for this circuit with an input 
at 300 baud is given in Fig. 14a. When the lock range was increased to 
±300 Hz about the carrier, the eye pattern improved, as shown in Fig. 
14b. This effect results from the increased bandwidth, / e . 

Start-stop distortion, as measured with a Bell System 911 A data test 
set for bit rates of 150 to 300 baud, was 2 and 5 percent, respectively. 
Performance with additive gaussian noise is as shown in Fig. 15, 




M 

m 




I 

■i 

ii 



(b) 



Fig. 14 — First order PLL eye patterns for 300 baud random digital data. 
Lock range of (a) 400 Hz, (b) 600 Hz. Horizontal scale: 0.5 ms per cm. 
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Fig. 15 — Distortion performance of first order PLL with input signal de- 
graded by 3 kHz band-limited white noise. 
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(b) 

Fig. 16 — Second order PLL eye patterns for (a) 150 baud and (b) 300 baud 
random digital data. Horizontal scale: 1 ms per cm. 
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where start-stop distortion is plotted versus signal-to-noise ratio. The 
noise is 3 kHz band-limited white noise and the bit rate is 300 baud. 
This performance is similar to that of the receiver in Data Set 103E1. 
Start-stop distortion for a single frequency (2025 Hz) interference 
was measured and, as might be expected, single frequency interference 
is a strong function of frequency, the worst values being at odd multi- 
ples of the channel frequency. 

Laboratory tests indicated that the second order PLL, with char- 
acteristics described in Section VI, performed in accordance with 
theoretical expectation. The eye patterns for an input of a modulated 
random data signal at 150 and 300 baud are shown in Fig. 16. It was 
found that noise performance of the second order system was not 
much better than that of the first order system. This probably resulted 
from the counter length falling just short of the value indicated in 
equation 78, thereby allowing noise perturbations to randomly shift 
operation between the stable and parasitic modes of operation. 

VIII. CONCLUSIONS 

A synthesis procedure for the realization of an nth order digital 
phase-locked loop has been described. Such systems find application 
in FM demodulators, filters, and in extremely stable locked oscillators. 
Various loop properties have been theoretically derived and experi- 
mental performance has been found to be consistent with these results. 
The advantages of such circuits for use in large multichannel data 
sets are: 

(i) Filtering property — (6n) dB per octive. 

(ii) Small size — completely integrable using one or more mono- 
lithic chips. 

(Hi) Requires no adjustment — permitting lower manufacture and 
repair costs. 

(iv) Excellent stability and reliability — the PLL circuit either 
works or does not, since it is completely digital. Stability of the en- 
tire system is dependent upon clock stability which may be as good 
as required. 

(v) Multichannel economy — accurate clocks can be used to drive 
many channel circuits. 

In addition, this circuit has two inherent advantages over other 
types of phase-locked loops. First, it requires no low-pass filter gen- 
erally found between the phase comparator (multiplier) and voltage- 
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controlled oscillator in conventional phase-locked loops. Second, it 
includes, in effect, an ideal voltage-controlled oscillator, the frequency 
of which is linearly related to voltage. 
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